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Abstract 

We consider Hermite and Laguerre /3-ensembles of large N x N random matrices. For 
all /? even, corrections to the limiting global density are obtained, and the limiting density 
at the soft edge is evaluated. We use the saddle point method on multidimensional inte- 
gral representations of the density which are based on special realizations of the generalized 
(multivariate) classical orthogonal polynomials. The corrections to the bulk density are os- 
cillatory terms that depends on (3. At the edges, the density can be expressed as a multiple 
integral of the Konstevich type which constitutes a /3-dcformation of the Airy function. This 
allows us to obtain the main contribution to the soft edge density when the spectral param- 
eter tends to ±oo. 
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1 Introduction 

We deal with two families of N x N random matrices: the Hermite and Laguerre /3-ensembles 
(for a review see (§]). These ensembles possess an eigenvalue joint probability density function 
(p.d.f.) of the form 



-PW(x) 



x = (xi, . . .,x N ) G I 



N 



(1) 



where (3 is real and positive. The support I of the eigenvalues in the Hermite and Laguerre 
cases are respectively (—00,00) and (0, 00). The ensembles' names come from the fact that their 
p.d.f. generalize the weight functions related to the Hermite and Laguerre polynomials; that is, 
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1 N 

E^- E ln 
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Hermite, 
Laguerre, 



(2) 



where a is a real and nonnegative parameter. The normalization constants can be computed 
with the help of the Selberg integrals: 



Z 



N 



N 



TT r(i+j/3/2) 

P' N =9^1} r(l + /3/2) 

3 — 2 



r(l+j/?/2)r(l + (a + i-l)/3/2) 

r(i + /3/2) 



Hermite, 



Laguerre, 



(3) 



where g , N = {2^/2 r N{x/2+p{N-i)/A) and = (2//3) iV(a^/ 2+ i +/ 3(Ar-i)/2) _ 

For special values of the Dyson index f3, we recover classical random matrix ensembles (see 
e.g. Indeed, the = 1,2, and 4 Hermite ensembles are respectively equivalent to the 
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Gaussian orthogonal, unitary, and symplectic ensembles. The Laguerre ensembles are similarly 
related to the real, complex and quaternionic Wishart matrices. Recently, Dumitriu and Edel- 
man pj have constructed explicit random matrices associated to the Hermite and Laguerre p.d.f. 
given in Eq.Q. A generic random N x N matrix belonging to the Hermite /3-ensemble can be 
written as a tridiagonal symmetric matrix: 



1 



/N[0,1] X(N-l)/3 

X(N-i)p N[0, 1] 

X(N-2)P 



V 



X(N-2)P 

N[0,1] 



X2/3 



X(N-3)/3 
N[0,1] 

Xf3 



Xp 

N[0,1] 



/ 



This means that the N diagonal elements and the N — 1 subdiagonal elements are mutually 
independent; the diagonal elements are normally distributed (with mean zero and variance 1) 
while the off-diagonal have a chi distribution. Recall that the densities associated to N[//, a] 
and Xk are respectively (27r<7 2 ) _1//2 e~ (a:_At < ) ^ 2<T > and 2x h ~ 1 e~ x /T(k/2), where in the latter 
case x > 0. Any N x N matrix Lp of the Laguerre /3-ensemble also has a tridiagonal form: 
for some N x N matrix 



Br 



1 



(xpp 



\ 



X{N-1)P 

X(p-i)p 



\ 



X(N-2)P 



P 



N + l--. 



X{P-N+i)p Xp J 



In this article, we compute the density for large but finite random matrices of the Hermite and 
Laguerre /3-ensembles. The density, or the marginal eigenvalue probability density, is defined as 
follows: 



PN,p(x) :z 



N 
Zn 



Pn,p(xi, ■ ■ .,x N )dx\ ■ --cIxn. 



(4) 



The quantity A r_1 pN,p(x)dx represents the probability to have an eigenvalue in the interval 
[x,x + dx). The density has two simple physical interpretations. 

First, we remark that the Hermite p.d.f. is equivalent to the Boltzmann factor of a log- 
potential Coulomb gas with particles of charge unity confined to the interval (—y2N, y2N) 
with neutralizing background charge density — (y2N /it)\Jl — x 2 /2N. From this point of view, 
Zn (divided by AM) is simply the canonical partition function at inverse temperature (5 and 
PN,p( x )dx gives the number of charges present in the interval [x,x + dx]. This analogy allows 
one to predict the global density : 



lim 

N— >oo 



PN,p( 




(5) 



This result is known as the Wigner semicircle law. For a finite matrix, we expect that the 
scaled density is of order one in the interval (—y/2N, s/2N), the 'bulk region' of the mechanical 
problem, while it decreases rapidly around ±v / 2A r , called the 'soft edges'. A similar log-gas 
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construction is possible for the Laguerre case. One expects the c = 1 Marcenko-Pastur law [181 : 



lim 4p N p(4Nx) = pup(x) 



r 2 



7T 



< x < 1, 
x > 1. 



(6) 



We see that, in the Laguerre case, the 'bulk' is (0, 4JV) while the 'soft edge' is the point 4iV. 
The origin is referred as the 'hard edge' of the support because the eigenvalues are constrained 
to be positive. The predictions given in Eqs (JSJ) and © have been confirmed in |UE|. The 
asymptotic analysis used in these references constitutes the starting point for the study of the 
higher expansions to be undertaken in the present work. 

Second, their is a deep connection between the /3-ensembles and some integrable quantum me- 
chanical iV-body problems on the line, known as the Calogero-Moser-Sutherland (CMS) models 
(a good reference is |21j). The Hermite p.d.f. is in fact the ground state wave functions squared 
of the (rational) Ajy-i CMS model, whose Hamiltonian is 



d 2 



N 

^dx- 



+ 



t 
4 



N 
i=l 



W - 2 ) 



l<i<j<N 



for Xj G 



oo,oo). The Laguerre p.d.f. is the ground state squared of the Hamiltonian of the 
Bn CMS model, which can be expressed as follows: 



N 



2 yi2x t — + — 

1=1 



1 



N 



8=1 



+ LX2[af){af3-2)-+f} 2 x i ) +/3(/3-2) ^ 



Xi + x 



3 



l<i<j<N 



(Xi 



where Xj £ (0,oo). It has been shown in P (see also [Hj) that the eigenfunctions of the 
conjugated Schrodinger operators e@ w ^ 2 e ~P w / 2 an d e@ w ^ 2 e~@ w ^ 2 are respectively the 
generalized (or multivariate) Hermite and Laguerre polynomials, previously introduced by Las- 
salle in [171 In the context of CMS models, the global density can be seen as the ground 



E 



j=l,...,N 



5(x — xj), also known as the 



state expectation value of the density operator p(x) 
one-point function. 

The relation between the CMS models and the generalized classical orthogonal polynomials 
furnishes, when (3 is an even integer, new integral representations of the global density that 
suits perfectly for asymptotic analysis. Let us be more explicit. The definition of the density 
given in (jlj) contains iV integrals; considering N large does not simplify the calculation. On the 
other hand, it has been noticed in ^[7] that the density is a particular Hermite (or Laguerre) 
polynomial, characterized by a partition A = {{N — If) and evaluated at x\ = . . . = xp = x 
(see below). Using the work of Kaneko ^1] and Yan [22j . one then can shows that the density 
is proportional to the following /^-dimensional integral: 



Rn,p(x) ■= / duie 



du 
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i<i<fc</3 



U.j 



Uk 



4/(3 



(7) 
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for a particular contour C and function f(u,x). 

In the following sections, we apply the steepest descent method [2UJ to integrals of the 
type • We obtain expressions for the density in the bulk and at the soft edge that are valid for 
every (3 S 2N. Of course, these results generalize many known result obtained for (3 = 2 and 4. 
We mention in particular two recent publications in which asymptotic corrections to the global 
density have been obtained: 1) Kalisch and Braak ^3] for some ensembles, including the Gaus- 
sian unitary and symplectic ensembles (work based on the supersymmetric method); 2) Garoni, 
Frankel and Forrester jllj for the Laguerre and Gaussian unitary ensembles (calculations using 
the theory of orthogonal polynomials). Also, the preprint [It), of Forrester, Frankel and Garoni 
addresses the Laguerre and Gaussian ensembles with orthogonal and symplectic symmetry. All 
studies show that these approximate expressions of the global density are very accurate, even for 
N = 10, say (for instance, see Fig. 1 and Fig. 2 in JJ)- We finally point out that an asymptotic 
formula for the density in the Hermite /3-ensemble has been considered in a different context: 
Johansson ^2j has studied a smoothed (macroscopic) density and has derived corrections of 
order 1/N to Eq. However, contrary to the asymptotic formula obtained here, the large iV 
expansion given in does not contain oscillatory (microscopic) terms. 

The article is organized as follows. In Section 2, we review the exact expressions of the 
densities in terms of the generalized Hermite and Laguerre polynomials. In Section 3, we derive 
the first oscillatory corrections to the global densities © and These approximations are 
also compared to the exact densities given in Section 2. The asymptotic densities evaluated 
about the soft edges of the spectrum are obtained in Section 4; they are expressed in terms of 
Kontsevich type integrals. The behavior of the latter when the spectral parameter is large is 
studied in Section 4. In the last section, we finally summarize the principal results and discuss 
the generalization of some of our results to general (3. 

2 Exact expressions of the density 

As previously mentioned, the density in the Hermite and Laguerre ensembles can be written as 
particular generalized Hermite and Laguerre polynomials pQ. These polynomials are symmetric, 
so we can write them as a linear combination of monomial symmetric functions 

m\(x\, . . . , xn) '■= x^ 1 ■ ■ ■ Xpf + distinct permutations, 

where A = (Ai,...,Ajv) is a partition of weight |A| = J2iLi^i- It * s convenient to introduce 
another basis of the algebra of symmetric polynomials, namely, the (monic) Jack polynomials 
jjf' . They constitute the only basis such that 

j{ a) (xi, . . .,x N ) = m A (xi, . . .,x N ) + ^2 a ^i a ) "vO^i, ■ ■ -,x N ) (triangularity) 

D2 J\{ x li ■ ■ ■ ,%n) = J\ . . . ,xn) (eigenfunction) 

for some eigenvalue e\(a). In the last equations, \i < A means that X^j=i Mi — Sj=i f° r an ^ 
when = |A| but \i 7^ A, while is a particular differential operator that can be defined via 

_(a) v^fcd 2 1 ( k 9 k 9 \ 

k ^ 0s? a , Xi - xj V dxi 3 dxj J 
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The generalized Hermite polynomials, denoted by H\(xi, . . . a), are the only symmetric 
polynomials obeying to 

H\(xi, . . . ,x N ;a) = J x \xi, ■ ■ ■ ,x N ) + ^ b Xfl (a, N) Jjf^ (x 1 , 

\[i\ = \\\-2n 
n=l,2,...,L|A/2[J 



(D { a) -2E 1 )H x (x 1 , 



,x N ;a) 



-2\X\H x (xi,. . . ,x N ;a), 



where 



N 



d_ 

dxi 



Let us point out that 



D, 



(a) 




2El = -leW 2 e-W 2 + est, 13=-, 
p a 



where is the CMS Hamiltonian defined in Section 1. On can show that 
H\(xi,...,x N ;a) = exp f~D^\ j[ a) 

Similarly to the Hermite case, the generalized Laguerre polynomials, written L^(xi, 
are the unique symmetric polynomials satisfying 

L u x (xi,...,x N ;a) = j[ a \x 1 ,...,x N ) + ^ c Xli (a, u, N) Jj^ (xi, . . . , x N ), 



\li\=\\\-n 
n=l,2,...,|A| 



(Z)j a) - Ei + + 1)E ) L v x {x\, ...,x N ;a) = -|A| L x (x 1 , ...,x N ; a). 



(8) 

,x N ;a), 



e= 2 -, 

a 



pa- 1 



The latter eigenvalue problem is related to a CMS model: 

D[ a) -E 1 + {u+ l)E = -^,e mt2 e~ m l 2 + est, 

Ip 

The following formula furnishes a way to compute the generalized Laguerre polynomials: 



L x (xi, ...,x N ;a)= exp ( -D^ - (u + l)E \ J ( x "'{xi, . . .,x N ). 



(9) 



When (3 is an even integer, the density in the Hermite and Laguerre ensembles can be respec- 
tively written as a particular Hermite and Laguerre polynomial; explicitly, 

G 

' < " -.11,, \ ,l.n r r. J/211 



PN,p{x) 



WV-l p -/3x 2 /2 j 



Hermite, 



N 



W a ,0, N -l a0/2 0x/2 fa-l+2/0 , XR .B/2)\ 



Laguerre, 



(10) 



where we have used the convention (n k ) = (n, . . . , n). Eqs Q and (jHJ), together with the fact 
that j£*\xi, . . . , Xk) = X™ ■ ■ ■ x^ when A = (n k ), readily imply the following exact expressions 
of the density: 



G 



P,N-1 -/3x 2 /2 



P(N-l)/2 



PnA x ) = N ~PT e 

G0,N 
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n=0 



4 n n\ 



o J x i 



(11) 
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in the Hermite case, and 



, (12) 

- X\=...=Xji=X 

in the Laguerre case. 

We will use the two latter formulas to compare the exact and asymptotic expressions of 
the density. Note however that Eqs Qllj) and Q12|) are computable only when both iV and 
(3 are small (i.e., for small partitions in Eq. There exist other methods that allow to 

calculate the multivariate classical polynomials (see for instance (Sj), but they suffer from the 
same restrictions. As shown in the next section, the asymptotic expressions provide a more 
tractable way to determine the density when TV is large. 



PnA x ) 



.w. 



a,f3,N-l ^g0/2 c -0x/2 



N ~^r — x 

0(N-l) 
n=0 



(-iy 



nl 



Df 2) + (a-2//3)E ) x{ 



JV-1 



3 Density in the bulk 

In this section, we obtain oscillatory corrections to the global densities ((HJ) and ©. This is 
achieved by deforming the contours of integration C in ((7J) in such a way that they pass through 
the saddle points of the function f(u,x). In both the Hermite and the Laguerre cases, the 
function f(u,x) has two simple saddle points in the complex u-plane, called u + and U-. All 
oscillatory terms can be seen as combinatorial corrections: the global density is recovered when 
(3/2 variables go through u + while the remaining (3/2 variables go through u_; the dominant 
oscillatory term comes from the integration of (3/2 + 1 variables through u + and (3/2 — 1 variables 
through U- and conversely; the second oscillatory term comes the integration of (3/2 + 2 variables 
through u + and (3/2 — 2 variables through u_ and conversely; and so on. 

3.1 Hermite case 

Before considering explicitly the density in the Hermite ensemble, we prove two technical results 
associated to the asymptotics of the integral (J7J). In the following lines, we suppose that e N ^ u ' x ' 
is analytic everywhere in the finite complex n-plane, except possibly at a pole depending on x, 
and that C is the real interval (— oo, oo). 

Figure 1: New contours {&., . . . ,C„} in the complex itj-planc. 

Uj + i 




The method of steepest descent requires that the integrand of Q should be analytic. This 
means in particular that the absolute values must be removed. Such an operation is realized in 
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the following lemma; it is possible when the line integration C of the variable Uj is deformed into 
an appropriate complex path Cj. Acceptable contours are given in Fig. ^ Other appropriate 
contours are obtained by making a reflection of the picture with respect to the real axis. Note 
that the dashed lines stand for (movable) branch cuts. We stress that Uj's contour starts at — oo 
and ends at the complex variable Uj-\. Only the path of u\ (the last variable to be integrated) 
ends on the real axis. 

Lemma 1. Let {Cj} be a set of non-intersecting contours such thatC\ is a simple contour going 
from — oo to oo and such that Cj goes from —oo to Uj—x for all j = 2, . . . ,n (see Fig. QJ). Then 



/OO fOQ n 
ah-/ ii' vr: "-" n 
-oo J —oo ■ i 1 ✓ ■ >l. 



i=l ^<3<k<n 

n 



[ dm- f du n f[e N ^ II ( Uj - u k fW, 

JCn i=l l<j<k<n 



where —ir < argUj < ir and where arg (n, — Uj) A /P = when Ui,Uj € R but Ui > Uj. 1 

Proof. Using the invariance of the lefthand side of the previous equation under any permutation, 
we immediately see that it is equivalent to the following ordered integrals: 

dui e Nf(u uX ) / ^ e Nf M . / ^ e Nf(u n , X ) -Q {u . _ Uk) A/P . 

-oo J-oo J-oo l<j<k<n 

These integrals only contain analytic functions, so we can use Cauchy's Theorem. This implies 
that the contour of ui can be deformed into any simple curve starting at — oo and stopping at 
oo. For the remaining variables, any nonintersecting contour of integration in the complex plane 
which starts at — oo, doesn't cross the branch cuts coming coming from the multivaluedness of 
the integrand, and complies with the ordering of the variables can be chosen. □ 

We are now in position to analyse ® when iV is large. Recall that the basic idea of the 
steepest descent method is to choose a path for which the decrease of /(«, x) in maximum. In 
particular, this means that the contour must pass through the saddle points. In the bulk case, 
/ has two simple saddle points u±; that is, 

|-/(n,x)| u± =0, ^- 2 f(u,x)\ u± =Re^, R>0. (13) 

The directions of steepest descent at these points, denoted 6±, are such that cos(2#± + <fi±) = — 1 
and sin(20± + <f>±) = 0, so 

0± = 71 ~ (modTr), -it<6±<-k. (14) 



1 When n 7^ (3 and (3 £ N in Lemma0 other integral representations in which all variables go from —00 to 00 
are possible 0- 
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Proposition 2. Let /(it, x) be a function that satisfies Eqs (|13|) and (|14|). Let aZso /-t 
f(u±,x). Suppose moreover that the saddle points are such that K(it_) < 5R(n + ). Then, 







(u+ — li_)^ 



— 1 



,3-1 



x e pN(f++f-)/2j(p-i){e++e-) 



rN,p( x ) + 



A/ 



where 



r N)f) (x) = 1 + 2 



fe=i 
e 



n T(l + 2j/f3) 

jlr(i + 2(j-k)/p) 



i2k 2 (e++8-)/f3 



(U+-U-) 



4k' 2 



and 



n,0 



du\ ■ ■ ■ 



l<j<fc<n 



i=l 



4//3 = _J 

2™( 



n/2 



n-l)//9 11 r(l 



j=2 



r (i + 2j-//?) 

(1 + 2/0) ■ 



(15) 



Proof. We first apply Lemma U to the expression (jJJ). Then, the contours Cj are deformed into 
steepest descent contours Sj = <S~ U 6+ passing through the saddle points u± . Close to these 
points, the contours are parametrized as follows: 



Uj = u± + tje 



-it < argnj < 7r, on Sj , 



where the angles of steepest descent are given ()14|) and where tj G (— r, r) for some r > 0. 
Moreover, we impose ij > ij for i < j in order to guarantee K(itj) > 9?(ity) for i < j. Setting 



NR 



Vi = 
we obtain 

Vl £ (-oo, oo) and x) = Nf ± - y) + 0(l/ViV) 

as N — > oo. When both itj and it^ are close to the same saddle point u±, 



i4fl±//J te-!/*) 4//J - 



When itj is on while it& is on <S • , we have 
(uj - u k flP = (u+ - uJflP + 0(l/\/iV) . 
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We now return to the expression of R^^{x) by considering the steepest descent paths: 



RnA x ) = P ] - yZ \ dm--- du n 



x du 



n+l 



Js d j=l l<j<k<0 J 



In terms of the new variables introduced above, the righthand side of the previous equation 
becomes 



13 



n=0 



yi 



dyi / dy- 



2 • 



Vn-l 



dy n 



+i 



Vn+l 



dy 



n+2 ' ' ' 



yp-i 



dy n 11 e-"i ( 1 + O 



l<p<g<?i 



fc=n+l 



* n 1+0 uw i n (*-^) 4/ " 



where 



5„ = (u+ - u_) 



4n(/3-n)//3 



2 \ 30/2-l-2n(p-n)/P 

Wr) 

x e JVn/++iV( / 9-n)/_ e ie + (n+2n(n-l)/^) e ie_(/3-n+2(^-n)(^-n-l)//3)_ 



The Gaussian terms are in fact ordered versions of the functions T n g, introduced in Eq . <|15(1 . 
We again use Lemma ^ and get 



R »d*) = £ (f) 

n=0 V y 



Sn^n,p^p-n,P + O ( — 



iV 



where 



/3/2 



/3 V /3 



1 r / g/2+fc,/3r / 3/ 2 -fc,/3 



S/3/2+k + S/3/2-k 



S. 



0/2 



(16) 



Note that the order of the 'analytic corrections' (i.e., coming from the expansion of /(%■, x) 
as a polynomial in yj of degree superior than 2) is now 1/JV rather than 1/y/N. This can 
be explained by using the theory of the generalized Hermite polynomials T: only symmetric 
polynomials p(yi, ■ ■ ■ , y n ) of degree even may have a nonzero contribution to 



/oo poo n 2 
dyi--- dy n ~[[e~ y i ] J - y fc | 4//3 p(?/i, . . . , y n ). 
-OO J— OO ■ -| T^.^l ■ 



i=l l<j'<fc<n 
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In our case, the 0(1/ yN) terms are symmetric polynomials in yj of degree one and three, so 
they don't contribute in the expression of Rjya(x). 

Note also that (S^/2+k + Sp/2-k)/Sp/2 * s of order 1/N 2k2 ^ . Hence, the 'combinatorial terms' 
with k > Lv /3/2J are smaller than the 'analytic corrections' of order 1/N; thus we must truncate 
Eq. l|Trj|) as follows: 

L\/^2J / a X / a s -It- n ro , c 



This expression contains only 'combinatorial corrections'. 
Finally, one readily shows that 

P \f P V 1 r /3/2+fc ^r /3/2 _ fci/3 _ i * r(i + gj/g) 

v /?/2 + A,,/ V/?/2y' r mp T mf3 ~ f = [ T(l + 2(j - fc)/0) 

and the proof is complete. □ 

The explicit link between the Hermite density and the integral of the type (j7J) has been 
obtained in pQ. It reads: 

p N AVWx) = l --^^(2Nf N l^e-^ 2 R Nfi (x) (17) 

if 

f(u, x) = -2u 2 + ln(hi + x) - — ln(iu + x). (18) 
Up to additive terms of order 1/N, the latter function has two saddle points 
u± = — (ix ± \Jl — x 2 ) . 

We see that 3ft (u + ) > 3ft only if — 1 < x < 1. According to the notation used in Eqs (fT3*|) 
and (|14|). we have 



1 iV-1, 2 /iV-l / ~\ 

f + = — — In 2 + x db i — — — arccos x — xy 1 — ar 

J 2 N \ N J 



and 

RfM± = Sy/l _ x 2 e i(^-arc S mx) ; 

where we have made use of 



arcsinx = — im(ix + yl — x 2 ) = tt/2 — arccos x. 

Note that the inverse trigonometric functions are defined on their principal branch; that is, 
arcsinx : [—1,1] — > [—tt/2, tt/2] and arccosx : [—1,1] — ► [vr,0]. Thus the angles of steepest 
descent are 

0± = ^f- arcsinx . 
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Stirling's approximation, 

T{y + z) = V2^e~ z z y+Z - 1/2 



1 + 



when z — > oo, 



immediately implies 

G^jV-l _ 2 W 2 - 1 /2 /3 /3iV/ 2 -/3 e /3iV/2 



G 



0,N 



nN pN/2+l/2 



1 + 



N 



The substitution of the above results in Proposition |2 provide the sought asymptotic corrections 
to the global density. 

Corollary 3. Let — 1 < x < 1 and let Pw{%) denote the (cumulative) probability distribution 
associated to the semicircle law given in Eq. i.e., 



f x x 1 
Pw(x) = / pw(t)dt = 1 H — Pw( x ) arccosx . 

J-l 2 7T 



(19) 



Then 



2 PN,p(\/2Nx) = pw(x)r N ^(x) + O ( — 



where 



1+2 ^ 



^ (7r 3 pw(^) 3 iV) 2fc2 /^ 



/or (p(x,/3) = yl — — J arcsinx. 

Let us consider only the very first correction to the global density: 



(20) 



l W (V2iVz) = pw(x) + O (1) + O 

2 T(l + 2//3) 1 



vr (vr / o w (x)) 6//3 - 1 iV 2 // 3 



cos(27riVPw(^) + yj(x,/3)) . (21) 



Up to a factor of order 1/N, the dominant oscillatory terms in the Gaussian unitary and sym- 
plectic ensembles are thus 



p N>j3 (V2Nx) - pw{x) 



cos(2irNP w (x)) , 



13 = 2. 



{ vrpw(x)V 2 iVV2 COS ( 27riVPw W + \ a^csin x), (3 = 4, 
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respectively. A direct computation shows that the non-oscillatory 0(1/N) term is exactly zero 
when (3 = 2. This implies that our result reproduces the Gaussian global densities previously 
obtained in |11| I13j. for the unitary case, and in |13j . for the symplectic case. 2 

The asymptotic expansion of the density for (3 = 6 is numerically compared with the exact 
one in Fig. [21 This picture shows that, even for a small iV's, Eq. (J20J) furnishes a qualitatively 
good approximation of the density in the bulk. Fig. 01 illustrates the behavior of pN,p(V2Nx) 
when the Dyson index varies. 

Figure 2: Comparison of the exact density 1111 . shown as a solid line, and the asymptotic density 1211 . shown 
as a dashed line, in the Hermite /3-ensemble for N = 7 and /3 — 6. 



3.2 Laguerre case 

The method used to evaluate the asymptotic behavior the the Laguerre density is almost the 
same as the one used in the Hermite case. All relevant differences originate from the correct 
contour that we must choose in Eq. 7 j: C starts at the point u = 1, turns around zero in the 
counterclockwise direction, and comes back to u = 1. In the following paragraphes, we briefly 
obtain the Laguerre version of LemmaHand Proposition^! We suppose that e N ^ u,x ' is analytic 
everywhere, except maybe on the interval [0, 1]. 

Lemma 4. Let {Cj} be a set of nonintersecting counterclockwise contours around the origin, all 
starting at Uj = 1, such that < arg(n n ) < . . . < arg(ui) < 2ir. Then 




x 




n 




Uj - U k 



A/13 



1=1 l<j<k<n 



n 






2 Note however the presence of a misprint in [T3] for the symplectic case. 
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Figure 3: Asymptotic density £fll in the Hermite /3-ensemble for N = 8 and (3 = 2,6, 10. 
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x 



2(n — 1) 

where f(u,x) = f(u,x) -— — lnu. 

Proof. We first set \uA = 1, i.e., Uj = e l6,J . C is such that 9j goes from to 2ir. The integral 
is completely symmetric so that we have n\ possible arrangements of the type < 0^ < . . . < 
@i„ < 27T. We choose < 9 n < . . . < 6\ < 2tt. In that case, 



n K-^r= n ( 2sin ^) 



4//9 



The righthandside of the previous equation can be written as 

I5 _ i _L w n n (tl( _„. )W . 

i=l i<i<fc<i 
We have proved that 



id U1 --- <fdu n f[e N ^ u ^ n 

c ^ c i=i i<j<k<n 

/ K|=1 d Ul ...du n \{e N ^ [J ( Uj - unflP 



.]Wn-l)//3 



0<arg(n n )<...<arg(ui)<27r 1=1 l<j<fc<n 



The integrand is analytic everywhere but possibly on the segment [0, 1]. Therefore, we can apply 
Cauchy's theorem and deform the paths on the unit circle into any counterclockwise contours 
Ci around zero and starting at u\ = 1 as long as the ordering < arg(u n ) < . . . < arg(ui) < 2ir 
is satisfied. □ 

Proposition 5. Let f(u, x) = f(u, x) — (2— 2/(5)N~ l lnu, where f(u,x) is the function appearing 
in the definition of Rn,/3(x), given in Eq. (J7J), and satisfying Eqs (fT3)) and (|14|). Let also f± = 



3.2 Laguerre case 
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f(u±,x). Suppose moreover that the saddle points are such that < arg(u + ) < arg(u_) < 2n. 
Then, 



^W = -G/ 2 )(%v) 2 K-«-)^(^) 



p-1 



x e /3AT(/ + +/_)/2 e i(/3-l)(0 + +0_) 



r N>p (x) + 0(l/iV) 



where 



r N ^(x) = l + 2 



k=l 



Urn 



r(i + 2j-//3) 



vi =1 



+ 2(j - fc)//3) 



i2fc 2 (e++e_)//3 



cos - ikN{f+ - /_) + - 0_)(3 - 2//3) 



{u + -u.Y k2 /P{NR) 2k2 /P 

Proof. We first use Lemma 01 The remaining steps are similar to those of Proposition [2j □ 

In reference [7j, the density of the Laguerre /3-ensemble has been written in terms of gener- 
alized hypergeometric functions |14j : 



PN,/3 



(x) = N Wa+ ^> N - 1 x a V 2 e-W 2 iF^i-N 



a,/3,N 



+ l;a + 2;ti,...,t /9 )| 



There exist integral representations of the generalized hypergeometric functions [21] ■ For our 
purpose, the appropriate integral formula can be found in Chapter 11 of [S]; one easily shows 
that 



lFl A+l + (n — l)/a;t u . . .,t n )\ h 

2nB ^ 

"_ ,p da ! . ■ ■ __ f| 



M n (A5,l/a)2vri// ni ' 2;ri 



n 

l<fc<Z<n 



where C is as previously described and where 

T(l + A + B-C + jC)T(l+jC) 



M, 



a r(i + a - c + jC)r(i + b - c + jC)r(i + c) ' 

This implies that the Laguerre density in the bulk can be recast in an integral of the form Q : 

N W a+2 ,f3,N-l (4 Nx )a(3/2 e -2/3Nx 



p N ^{ANx) 



{2Tdf W a ^ N M p {a + 2/(3 - 1, N - 1, 2/(3) 



Rn,p{x) 



22 



provided that, in Eq. (J7J), C is a counterclockwise closed path around the origin and starting at 
u = 1, and 



1 / 2\ 

/(u, a?) = 4xu — Inn + ln(l — u) + — \ a — 2 + — ) ln(l — u). 



(23) 
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Neglecting factors of order 1/N, we see that this function has two simple saddle points, 



u± = - l±iW--l , 



which satisfy arg(n + ) < arg(u_) only if < x < 1. This implies 



c) 



— 2 f(u,x)\ u± =Re^, R = Wx^--l, 



1 



L ^ 



hence, the directions of steepest descent are #_|_ = 37r/4 and = 7r/4. We also have 
a + 4//3-4 



/± = 2x 



A/ 



In ± 2i I x\ 1 — arccos \fx =F i-r? arccos V%- 



The Stirling approximation readily gives 



W a+ 2,/3,N-l _ ( P 



1+ap 



N af3/2_ 



r(i + /3/2) 



W a ,/3,iv V2 7 r(l + o/3/2)r(l + (a + l)/J/2) 

Moreover, by using the Gauss multiplication formula, which reads 

ra-1 

iQ r(« + j/n) = (2vr)( n - 1 )/ 2 n 1 / 2 - n2 r(nz), 
i=o 

one can show that 

P 



N 



1 + 0l N 



( r /3/2,/3) 



0/2J M p (a + [3/2-1, N -1,2/(3) 

1 fp\ ^~ a/3 r(l + ap/2)T(l + (1 + a)(3/2) 



7T 



1+01 



2; r(i + f3/2)N a P+ 2 -v 

By substituting the above equations in Proposition [SJ we obtain the following result. 



(24) 



Corollary 6. Let < x < 1 and let Pmp(x) denote the (cumulative) probability distribution 
associated to the density pmp defined in Eq. ©; i.e., 



f x 2 
Pmp(x) = / pMp(t)dt = 1 + xpmp(x) arccos x . 

J0 7T 



(25) 



Then 



4p Ntf 3(4:Nx) = pup{x)r Nj p{x) + O 



where 



r NJi (x) = 1 + 2 



-1) 



^ (2ir 3 x 2 pMp(x) s N) 2k2 /P 
k 



X 



n r(i + 2j-//3) 
l = lr(i + 2(j -*)//?) 



cos ( 2irkNPyip(x) + k(p(x, (3) ) , 
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Figure 4: Comparison of the exact density 1121 . shown as a solid line, and the asymptotic density 1271 . shown 
as a dashed line, in the Hermite /3-ensemble for N = 4, (3 = 6, and a = 0, 1. 




0.2 0-4 0.6 O.S 1 

x 



ip(x,(3) = fl — — ^ — — 2a arccos (26) 
Again, we look at the first correction to the global density: 
4p N> p (ANx) = pup (x) + O ( — J + O ' 



nj ym/f 3 

2T(1 + 2/13) 



- ^^U^ cm ^ NPupix)+vix ^ (27> 

Hence, neglecting the possible factor 0(1/A^), the first oscillatory corrections to the global 
density in the complex and quaternionionic Wishart ensembles are 

— 2 — cos (inNPywix) — 2a arccos \/~x) , f3 = 2, 



Ap N ^(ANx) - pmp{x) 



it 3 x 2 pmp(x) N 

-7- cos ( 27tA^Pmp(^) — 2a arccos J~x + 7r/4), 8 = 4. 

l2 1 /2 vra; p MP ( x ) 1 /2 7V i /2 v ; 



Contrary to the Gaussian unitary ensemble, the first correction in the complex Wishart ensemble 
has a non-null (when a / 0) correction of order 1/N which is not oscillatory Nevertheless, 
our oscillatory term is the same as the one given the latter reference. The = 4 case has been 
recently studied in JHJ; the dominant correction is purely oscillatory and is equal to the term 
given above. 

Fig. |U provides a numerical comparison between the asymptotic and the exact expressions of 
the density. Clearly, the asymptotic approximation is better for a = 0. Note that the oscillations 
shift to the right when a increases. In Fig. 03 we illustrate the effect of the variation of (3 on 
4p^ : p(ANx) for fixed N and a. 



4 Density at the soft-edge 
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Figure 5: Asymptotic density in the Laguerre /3-ensemble for N = 5, a = 0, and j3 = 4, 6, 8. 




X 



4 Density at the soft-edge 

We have seen that the steepest descent method can be applied to the scaled densities only if 
— 1 < x < 1 (Hermite case) or < x < 1 (Laguerre case). Indeed, when the spectral parameter 
x is outside these intervals, the contours of integration cannot be deformed into the steepest 
decent ones without transgressing the appropriate ordering of the variables of integration. A 
change of scaling is mandatory. The appropriate changes of variable at the soft edges have been 
obtained in [£]: the scaled densities pn,/3(V2Nx) (Hermite) and pn^(4Nx) (Laguerre) should 
be replaced by p Nt/3 (V2N + X/V2N 1 / 3 ) (Hermite) and Pn^N + 2(2A0 1/3 x) (Laguerre). 

Technically, these new scalings make the two simple saddle points coalesce and become a 
double saddle point (or saddle point of order two). Then, the multiple Gaussian integrals are 
replaced by multiple Airy integrals, or integrals of the Kontsevich type |1F 



1 pioo plOO 

KnA*) ■= "tA^ / dm- dvnH ]"[ \v k - Vl \^ . (28) 

{ 1 J - io ° J - lo ° j=l l<k<l<n 

We recall that the Airy function of a real variable x can be defined as follows: 

-i rioo 

Ai(x) = — / e v / 3 ~ xv dv, (29) 

27T1 7-ioo 

and as a consequence 

Ai (x) = xAi(x). 
One readily verifies that 



K 1>/3 (x) = -Ai(x), K 2j2 (x) = 2(Ai'(x) 2 - xAi(x) 2 ). 



4.1 Hermite case 
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It is worth mentioning that the function K n ^{x) has previously been studied in the context of 
the Gaussian unitary ensemble |5j- In particular, it has been shown that 3 



Kn,2( x ) = ~ n ^- det 
4.1 Hermite case 



d i+ i- 2 
dx i+ i~ 2 



Ai(x) 



i,j=l,...,n 



The next lemma generalizes a basic fact of the Airy function of a complex variable z; that is, 
the contour in Eq. (|29|) can be deformed so that 



Ai(*) 



1 

2vri 



e v 3 /3-zv dy 



Ao 



(30) 



where .4.0 is a simple path going from ooe l9a to ooe ldb , where —tt/2 < 6 a < —tt/6 and tt/6 < 
9 b <TT/2 (see in Fig.©. 

Figure 6: Possible contours in the integral representation of the Airy function. 
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■3tv 
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Lemma 7. Let Ao be the contour described above. Let also {Vj} denote a set of non-intersecting 
paths such that V\ = Ao and such that, for all j £ {2, . . . ,n}, Vj follows Aq but stop at Vj—\. 
Then 



IV. 



o v?/3-xvj 



(27Tl)-7 Vl J Vn f = \ xsk<l ^ 



n (vk-v^, 



(31) 



where —it < argVj < tt and where arg («j — Vj) A /^ = when both Qvi = = Qvj and Kw, > $lvj. 



3 Note however that the term n\ is missing in 
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Proof. We essentially proceed as in Lemma ^ Firstly, we write Eq. 1)28(1 as ordered integrals 
along the imaginary axis and remove the absolute values. Secondly, the analyticity of the 
integrand and the property 

lim [e" 3 /3-H =0 if -1<9<-^ *<9<*, or ^<0<^ (32) 

R^oo I \ v =Re^ 2 6 6 2 6 Q y ' 

are used to deform the contour of v\ into Aq. We finally complete the proof by exploiting 
the ordering of the variables, the choice of principal branch for (vi — Uj) 4///3 , and the Cauchy 
Theorem. □ 

Proposition 8. The integral Rn,/3 defined by Eqs (0 and ()18|) satisfies 
Proof. According to Lemma ^ we have that 

i 1 + ^m) = / ■••/** n ^ /k ' i+ ^ /2jv2/3) n k - ^ 



Considering x/2N 2 / 3 <C 1, 



i=l l<j<k</3 



f[^ 1 + W^)=9(u 1 x) + ——^-— N + 



2AT2/3; ^ ' > 2(m + x)Ar 2 /3 ^4/3 
where 

g(u, x) = — 2u 2 + ln(iu + x). 
This functions has a double saddle point uq = i/2: 



giu,x J 



du 



.2 



= = — -^g(u,x) 

u=u du 



d 3 

J)Z3S(U,X) 
u=u ua 



Re^o = l6e~ i7r/2 . 



The directions for which the decrease of / is maximum are determined by both conditions 
cos(3#o + 4>o) = — 1 an d sin(36*o + <po) = 0. Thus, the angles of steepest descent are 9o = 
— 57r/6, — vr/6, 7r/2. We choose the two former and make the following change of variables: 

vj = 2iN 1/3 ( Uj -u ). 

This implies that 

Nf 0* 1 + 2^) = f " Nln2 + ^ + ln2 + ^ " XVj + ° (jV^Ts) ' 

Let {Vj} denote the set of ordered and non-intersecting contours of steepest descent: Vi starts 
at ooe~ m / 3 , passes through the origin and stops at ooe 17r / 3 ; Vj follows Vi but stops at Vj-i, 
where j = 2, . . . , (3. We have proved that 



?!- 



e 



/3N/2+/3N 1 / 3 o 



L**"L**w-' n fe-^) 4/ "+ 0(^75) 



V /5 i=l l<i<fc</3 
Consequently, we can apply Lemma [7| and the proposition follows. □ 
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Let us go back to Eq. (|17[) and make the change of scaling xwl-f- x/(2iV 2 / 3 ): 
P7Vi/3 ( V2N + 



V2iVV3. 

' G m^{ K 2Nf N l 2 ^e^ N e-^ /3x e-^ 2 / { - ml '^R N ,p ( 1 + 



Proposition IHJand some manipulations directly imply the following. 

Corollary 9. The density in the Hermite (3-ensemble evaluated at soft edge is proportional to 
the an integral of the Kontsevich type: 

1 =p N0 ( V2N + 



(3/2 



i (±k\ p " r(i + /?/2) / l \ 

U y n? =2 r(i + 2/ / ?)- i r(i + 2,7/3) AM j + WW ' 



4.2 Laguerre case 

It is obvious from Eq. (|30|) and Fig. El that 



1 

2vri 



e v*/3-zv dv _ q ^ e - i27r / 3 Ai(e- i27r / 3 z)+Ai(z)+e i27r / 3 Ai(e i27r / 3 z) = 0. 



This result can be extended to the Kontsevich type integral as follows. 

Lemma 10. Let A-i,A± be the contours depicted in Fig. Let also {Vj} denote a set of 
non-intersecting paths such that Vi = A-\ U A\ and such that, for all j £ {2, . . . , n}, Vj follows 
Vi but stops at Vj-i. Then 

Kn^x) = (-l) n+1 7^k L dvi--- [ dv n n Yl {vk - Vl )W, (33) 

^ ' JVl JVn j=l l<k<l<n 

where —tt < argVj < ir and where arg (vi — Vj) A /^ = when both Qvi = = ^Vj and $lvi > 3lvj. 

Proof. By virtue of Eq. (|3*2*)) . it is possible to join A-i and A\. Thus, Cauchy's Theorem implies 
that Vi can be replaced —Ao- The constraint on the ordering of the variables {v2, ■ ■ ■ ,v n } is 
then considered and equivalence between Eq. (|31|) and Eq. (|33|) follows. □ 

Proposition 11. The integral Rn,/3 defined by Eqs (J7J) and ()23|) satisfies 

( X \ 2pN e f3{2NY^x / l \ 

r »j> y 1 + pNfjs ) = ^ [ r^mm^ K ^ Ax) + [w^ ) ■ 

Proof. We essentially follow the proof of Proposition |H1 By virtue of Lemma HJ we have that 

rn, p ( i + ^3 ) = I • • • 1 n n («,- - n fc) ^ ; 



i=l l<i<fc</3 
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where 

/ («> 1 + (^) = + (2^73 - + ^ (° " 2 + I) 111(1 " U) + ^ G " 2 ) ln(U) ' 

for 

g(u) = 4u + ln(l — u) — ln(u) . 

The latter function possesses a double saddle point uq = 1/2. One can check that 

d 3 g(u)/du 3 \ u=U0 = 32e i7r , 

so the steepest descent angles are 6*o = 0, 27r/3, Att/3. The contour of ui is chosen such that: (1) 
it approaches no by following the real axis in the negative direction; (2) it leaves the saddle point 
with an angle 2ir/3; (3) it turns around the origin in the positive direction; (4) it comes back to 
the uq with an argument of 4-7r/3; (5) it finally leaves this point and reaches the point u = 1 by 
following the real axis. The ordering of the variables around the origin implies moreover that 
m follows U\ but stops at Uj_i. When N is large, step (3) is irrelevant and the steepest descent 
contour brakes into two disjoints paths, namely (l)-(2) and (4)-(5). Now we set 

Vj = 2(2A0 1/3 (u -^), 
which means that 

Nf («,-, 1 + ^) = 2N + (2iV)V3 x + (4 - | - a) ln2 + - ^ + O (^) . 

When N — > 00, the contours of the variables {dj} that give the major contribution to the 
integral, denoted by {Vj}, behave as follows: Vi is the union of the path that begins at —00, 
passes close to the origin and ends at ooe~ 17r / 3 together with the path that starts at ooe 17r / 3 , goes 
near the origin and stops at — ooe -171 "/ 3 ; Vj follows Vi but stops at Vj-i, where j = 2, . . . ,{3. 
Therefore, 

e 2/3N e P(2Ny/ 3 x 



2AT)2/3 J ^-2^/3+4/3^/3-2/3 
X 



JVl Jv P i=i i<i<fc</3 



Lemma ITUl finally provides the sought for result. □ 
We apply the last proposition to the scaled expression of the Laguerre density given in (I22|) : 



PN,f3 (4N + 2(2N) 1 / 3 x 



(2vri)/3 Mf,(a + 2/p-l,N-l,2/P) N ' P \ + (2N) 2 /z J ' 

Minor manipulations and use of Stirling's approximation give the sought limiting soft edge 
density, which is identical to that obtained in Corollary 03 for the Hermite /3-ensemble. 
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Corollary 12. The density in the Laguerre (3-ensemble evaluated at the soft edge is proportional 
to the an integral of the Kontsevich type: 

2{2N) 1 ^p N ^ (AN + 2(2iV) 1/3 

(3/2 



-1 f^LY r(i + /3/2) / i \ 

2tt U J UU + 2//3)- X r(l + 2j/(3) M[X) + WW ' 



5 Asymptotics of the Kontsevich type integral 

Here we consider the leading order of Kp g(x) when x — > ±oo. This allows to match the soft 
edge density with the bulk density expanded about the edge. 

Proposition 13. When x is large and positive 

r e -fx 3 / 2 / 1 



(2ir)P x 3/3/4-i/2 ^3/9/4+1 

Proof. Following the discussion in the proof of Proposition |S1 we first change the contours in 
the Kontsevich like integral: 

^ > JV i JV P i=l l<j<k<0 

where {Vj} is such that V\ goes from ooe~ ie to ooe , passing through the point -^/x, and such 
that Vj goes from ooe~ ld to Vj-i for all j = 2, . . . , ft and 7r/6 < 9 < tt/2. We now set 

Thus, 

i=l l<j<k</3 

where Wj = e _17r / 2 Vj. By virtue of Lemma ^ we can write 



l^TTj X 7_oo J_oo . =1 l<i<fc</3 







-w*+0(*- 3 / 4 ) TT I,.,. _„„|4//3 



Note that the term 0(x -3 / 4 ) is odd in lOj. As explained in the proof of Proposition 1171 this 
implies that the actual correction to the integral is of order x _3//2 . We finally obtain the desired 
expression by comparing the last equation with Eq. (|15l) . □ 

Applying Proposition 1131 the next result gives the behavior of the density when the spectral 
parameter leaves the bulk. 
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Corollary 14. Let o~(x) denote the density evaluated at the soft edge: 



a(x) 



lim om a ( V2N H , X ) , (Hermite) 



— »oo 

Then, as x — > oo 



lim 2(2N) L/6 p N/ 3 4iV + 2(2iV) i/3 x , (Laguerre). 

k iV— >oo ' ' 



2vr (4/3)/ 3 / 2 a;3j9/4-l/2 ~ ^3/3/4+1 

When the density is evaluated at points inside the bulk but close to the edge, we should 
observe both decrease and oscillation (see Fig. 12103). This is confirmed in next paragraphs. 

Proposition 15. Let x = —\x\. When \x\ is large, 

u , , (I>/2,/?) 2 ( P \ rr-r ( 1 



vr/ 3 \[3/2j v 1 ' x ' p \x 5 / 2 
where 



kn ~ 1+2 V IT r(1 + 2j//?) cos {%x\^ -l k (l-?W 

k x3 - 1+2 ^ 26fc2//3 | X | 3fc2//3 ii r(1 + 2(j _ k)/(3) cos ^ 3 fi 2 k y 1 )) 



k=i 1 1 \j=i 

Proof. By rescaling the variables in Eq. I|28|). we get 

-1 /-ioo /■ioo J 



3=1 l<k<l<p 



where 



f( v ) = ^v 3 + v. 

The function / possesses two simple saddle points, namely, v± = e ±17r//2 . We have f± = f(v±) = 
±2i/3 and f"(v±) = 2e ±17r//2 ; whence the angle of steepest descent are 9± = it/ 2 + n/4. The 
remainder of the proof is a straightforward application of Proposition El □ 

Corollary 16. Let o~(x) be the quantity defined in Lemma\ll\ When x — > —oo, we have 



i n /Tl r(l + /?/2) (A. |3/2 7T / 2 



+ — ^ +0 



5/2 y ^| x |6//3-l/ 

The previous result can be obtained directly from the asymptotic density in the bulk of 
the Hermite (or Laguerre) /3-ensemble. Indeed, the change of variable x i— > 1 — |x|/(2iV 2 / 3 ) 
in Eq. (|20|) and the development of this expression for A r_1 / 3 |a;| -C 1 reclaims Corollary 1161 
However, it is impossible to derive Corollary ^] from Eq. Q2UJI by such an expansion. Note 
finally that Corollaries 1141161 imply that the density at the soft edge of the Laguerre /3-ensemble 
is independent of a when both N and \x\ are large. 
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6 Concluding remarks 

The aim of the article was to determine the large- N asymptotic expansion of the density in the 
Hermite and Laguerre /3-ensembles when (3 £ 2N. 

We have shown that the first correction to the global density is purely oscillatory when (3 > 2 
and is of order N 2 ^ . In the Hermite ensemble of NxN random matrices, the density contains N 
peaks; the greater is (3 and the higher are the oscillations. The influence of the Dyson parameter 
on the oscillations is the same in the Laguerre ensemble. However, the density in the latter 
ensemble contains N — 1 summits and a (delta) divergence at the origin. 

These results agree with the large- (3 asymptotic analysis realized recently in jlj. More pre- 
cisely, it has been proved that for (3 — ► oo, the density in the bulk of the Hermite ensemble can 
be written as a sum of N Gaussian distributions centered at the zeros of an Hermite polynomial 
of degree N (and similarly for the Laguerre case). These conclusions are, of course, coherent 
with the log-gas analogy presented in Section 1. Note that no constraints on (3 are imposed in 
Consequently, we may surmise that our asymptotic formulas (|21|) and (|27|) are valid for any 
real (3 greater that 2, though the general method to prove this is still missing. 

We have also shown that the density of the Hermite and Laguerre ensembles are both pro- 
portional to a Kontsevich like integral Kp t p{x) when evaluated about the edges of the spectrum. 
Although the exact densities of the Hermite and Laguerre ensembles are quite different, the 
asymptotic analysis of Kp t p{x) has revealed that they approach the same function in the soft 
edge scaling, thus verifying the expected universality. The Kontsevich like integral itself is a 
special function generalizing the Airy integral and, as such, is worthy for independent study. 
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